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o 

£\J ' There are various analytical approaches to the mean electromotive force £ = (u X b) crucial in mean— field electrodynamics, with u 
and b being velocity and magnetic field fluctuations. In most cases the traditional approach, restricted to the second-order correlation 
^ approximation, has been used. Its validity is only guaranteed for a range of conditions, which is narrow in view of many applications, e.g., 
O , m astrophysics. With the intention to have a wider range of applicability other approaches have been proposed which make use of the 
■ so-called r-approximation, reducing correlations of third order in u and b to such of second order. After explaining some basic features 
' ■ of the traditional approach a critical analysis of the approaches of that kind is given. It is shown that they lead in some cases to results 
1 which are in clear conflict with those of the traditional approach. It is argued that this indicates shortcomings of the T-approaches and 
^s^j ' poses serious restrictions to their applicability. These shortcomings do not result from the basic assumption of the t— approximation. 
' Instead, they seem to originate in some simplifications made in order to derive £ without really solving the equations governing u and 
b. A starting point for a new approach is described which avoids the conflict. 

> 

vo 

1 Introduction 

o ■ 

. In mean-field electrodynamics the mean electromotive force £ = (u x b) due to the fluctuations u and 
b of t he fluid velo city and the magnetic field plays a crucial role ( Krause and Radler 1971 . 1980l . MoffattI 
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Il978l . lRadler1l200nl V A central problem is the determination of £ for a given motion as a functional of the 



q , mean magnetic field. Various methods have been used for that. 

One approach, which we call "traditional appr oach" or "approach (i)" in the following, was established 
c/2 . together with mean-field electrodynamics at all ( Krause and Radler 1971 . 1980l ). Most of the calculations 



of £ have been done on this basis using the so-called second-order correlation approximation (SOCA) or, 
what means the same, first-order smoothing approximation (FOSA). This approximation in its original 
form, that is, applied in the case of a purely hydrodynamic background turbulence, ignores all higher than 



5-h ' second-order correlations in the fluctuations u of the velocity field (see section f3.ip . It can be justified 
i only in cases in which these fluctuations are not too large. The usually given simple sufficient conditions 
for its validity are in view of astrophysical applications rather narrow. Basically it is possible to proceed 
to higher-order approximations but this requires tremendous efforts and has been done so far only in a 
few simple cases (Nicklaus and Stix 1988, Carvalho 1992, 1994, Radler et al. 1997, see also section [5Tj) . A 
slight modification of the second-order correlation approximation applies also to the case of a magneto- 
hydrodynamic background turbulence (see section EPj) . 

In some more recent investigations (which are cited below) other approaches are used, which rely in 
a sense on the r-approximation of turbulence theory ( Qrszag 197Cll ) and are called "r-approaches" or 



"approaches (ii)" in the following. They go in so far beyond the second-order correlation approximation 
as they consider also higher-order correlations, which are then in the sense of a closure expressed by 
second-order ones. In approach (i) the relevant equations are simplified by a well-defined approximation 
and then really solved, and £ is calculated with these solutions. In the approaches (ii) a relation for £ is 
deduced from the original equations, but without really solving them. Instead, in order to get manageable 
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results, assumptions on the connection of some of the occurring quantities with £ are introduced. The final 
result for £ is to a large extent determined by these assumptions. The approaches (ii) cover from the very 
beginning also the case of a magnetohydro dynamic background turbulence. 

The main purpose of this paper is a critical analysis of the approaches (ii). Each step of approach (i) 
can be justified by the underlying induction equation or, in the case of a magnetohydrodynamic back- 
ground turbulence, induction equation and momentum balance. Clear, at least sufficient conditions for 
the applicability of the second-order approximation can be given. There is no doubt in the correctness 
of its results in the so defined range of applicability. It seems reasonable to assume, and we do so in this 
paper, that there is at least some overlap of the ranges of applicability of the approaches (i) and (ii). 
We have then to require that in thes e overlapping ranges the results of both approaches coinc i de. S imple 
versions of approach (ii) as used by IVainshtein and Kichatinovl ( 19831 ). Blackman and Field (|2002h and 
Brandenburg and Subramanian ( 2005bl ). called "simple r-approach" or "approach (ha)" in the following, 
deliver results which do not in all ca ses satisfy thi s requ irement. As we will show below the more s o phist i- 
cated version used in the papers bv lRadler etail <j2003h and bv iRoeachevskii and Kleeorinl (|200l 120041 ) . 
called "spectral r-approach" or "approach (hb)", does not satisfy this requirement, too. We have to con- 
clude that these approaches are not in full accordance with the basic equations mentioned. Therefore the 
results can not be taken for granted. We will propose a starting point for an alternative approach which 
avoids conflicts with approach (i). 

In section [2] we define the frame of our considerations and deliver the basic equations. In section [3] we 
recall the fundamentals of approach (i) and review some of its basic results. In section [3] we explain the 
approaches (iia) and (iib), derive a few results, restricting attention to the simple case of a non-rotating 
fluid, and pinpoint shortcomings of these approaches and deviations of the results from those of approach 
(i). In section [5] we explain a proposal for the alternative approach mentioned. Finally in Section [6] we 
summarize our findings. 



2 The mean field concept in magnetofluiddynamics 

Let us first define the frame of our considerations. We want to explain essential features of the various 
approaches to the mean electromotive force £ but do not strive to a high level of generality. 



2.1 Mean-field electrodynamics 

Let us consider the magnetic field B in an electrically conducting moving fluid. We assume that it is 
governed by the induction equation 

d t B = r ? V 2 B + V x (U x B) , V-B = 0. (1) 

Here U means the velocity of the fluid and r/ its magnetic diffusivity, assumed to be constant. Until further 
notice we consider the fluid motion, that is U, as given. 

We assume that the motion and therefore also the magnetic field shows irregular, e.g. turbulent, features. 
Any field F of this type is split into a mean field F and a "fluctuating" field /, that is, F = F+f. The mean 
field F is defined as an average of F. It is assumed that the averaging procedure satisfies the Reynolds 
averaging rules. Alternatively to the notation F we use also (F) in the following. 

Averaging the induction equation ([1]) we obtain the mean-field induction equation 

d t B = t]V 2 B + V x (U x B + £) , V • B = (2) 
with the mean electromotive force £ given by 



£ = u x b . 



(3) 



, 2008 11:45 Geophysical and Astrophysical Fluid Dynamics SOCAetcl5 



Mean-field electrodynamics: Critical analysis of various analytical approaches 3 

Elaboration of mean-field electrodynamics means studying the properties of this quantity. For this 
purpose we rely on the equation for the fluctuations b which follows from ([T]) and ([2]), 

8 t b = r]V 2 b + Vx(Uxb + uxB) + G, V • b = . (4) 

Here G stands for a term of second order in u and b, 

G = Vx(uxb)', (5) 



with (u x b)' = u x b — u x b. 

As can be concluded from ([3]) and @ the mean electromotive force £ depends, apart from initial and 
boundary conditions for b, on U, u and B. More precisely, if is a functional of these quantities in the 
sense that £ in a given point in space and time depends on U, u and B also in other space points and 
at past times. This functional is linear in B. For most of the applications it is reasonable to assume that 
it is in addition local in the sense that £ in a point in space and time depends only on U, u and B in a 
certain surroundings of this point. Likewise in many applications the variations of B in space and time 
are sufficiently weak so that the behavior of B inside the relevant surroundings of a given point can be 
represented by B and its first spatial derivatives in this point. Then we have 

£ l = £^ +a ij B J + b tjk -^-, (6) 

where the quantities £*> , djj and are functionals of U and u in the above sense but do not depend 
on B. 

Assume for a simple example that there is no mean motion, U = 0, and that u corresponds to a 
homogeneous isotropic, not necessarily mirror-symmetric turbulence. Since no isotropic vector exists, £^ 
must be equal to zero. Further a%j and have to be isotropic tensors, = a <5jj and b^, = (3 e^, with 
factors a and (3 being averaged quantities depending on u. Consequently we have 

£ = aB- (3V xB. (7) 

For another simple example we remain with U = and assume that u corresponds to an inhomogeneous 
turbulence in a rotating system, that is, under the influence of a Coriolis force. More precisely, u deviates 
from a homogeneous isotropic mirror-symmetric turbulence only by an inhomogeneity and therefore an 
anisotropy described by a vector g parallel to the intensity gradient, and by the anisotropy and the 
deviation from mirror-symmetry described by the angular velocity ft responsible for the Coriolis force. 
Then £^ needs no longer to be zero, and a%j and bijk are no longer purely isotropic tensors. For the sake 

of simplicity we nevertheless at first ignore £- ' , referring to the comments in sections El 13.1-H 13.1.51 13.2.11 
and 15.11 We further assume that the influences of g and ft on the turbulence are weak enough so that £ 
is linear in these quantities. Considering then the symmetry properties of the basic equations (see, e.g., 
Krause and Radler 1971, 1980) we have 

£ = -7g X B - ai (g • ft)B - a 2 (g • B)ft - a 3 (O • B)g 

-/5 V x B-Si (ft- V)B-a 2 V(n-B), (8) 

where the coefficients 7, ax, ... 62 are determined by u only. (The choice of the signs in (JSj) , which follows 
some other papers, is however unimportant at this place.) 

Given the structures of £ as in these examples the only remaining task is the calculation of coefficients 
like a, f3, etc. The approximations we want to discuss in this paper concern only these calculations. 
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2.2 Dynamical aspects 

Later we will relax the assumption that the fluid motion, that is U, is given. In preparation to this we 
provide here the relevant relations. Restricting our attention to an incompressible fluid we assume that U 
is governed by the momentum balance, 

g(d t U + (U- V)u) = -VP + gvV 2 U- 2gQ x U + -(V x B) x B + gF , V • U = . (9) 

Here g is the mass density, P the hydrodynamic (including centrifugal) pressure, v the kinematic viscosity, 
considered as constant, and F means an external force. A rotating frame of reference is assumed, with f2 
being the angular velocity that defines the Coriolis force. 
Taking the average of equations §§§ we obtain 

g (d t TJ + (U • V)U) = -VP+ giyV 2 TJ-2gfl xU + -(V x B) xB + g(F + T), V • U = 0, (10) 



with a mean ponderomotive force T given by 



1 



jr = -(u-V)u + — (V xb) xb. (11) 
fig 

In view of both £ and J- the equations for the fluctuations u are of interest. They can be concluded 
from (191) and (fTOj) . 

d t u + (U- V)u + (u- V)U 

= -- V(p+ -B -b) + vV 2 u- 2fl x u + — ((B • V)b + (b • V)B) + f + T, V-u = 0. (12) 
g li fig 

T stands for the terms of second order in u and b, 

T = -((u • V)u)' - -^V(b 2 )' + -((b • V)b)' , (13) 
Ifig fig 



where ((u • V)u)' = (u • V)u — ((u • V)u) etc. 

For many purposes it is neces; 
infinitely extended fluid we find 



For many purposes it is necessary to eliminate the pressure term | V(p + ^B • b) in (I12j) . Assuming an 



with 



d t u = i^V 2 u + R + L + U + f + f , V-u = 0, (14) 



R = 2(n.v)v, */>(*) = ±- f v ' x u( \P dV 

47T X - X'l 



L = — ((B • V) b + (b • V) B + VM) , M(x) = — / n " , 3 -f^-U- d 6 x 

fig 2vr |x-x'| 

u = - ( u.v)u-(u.v)u-vy, v(x) = i- / d lMI^l : y)/a*' d 3 x , (15) 

27r Joe, l x - x I 

f = f + V 5 , g{*) = ±f ^-^dV 
4vr |x - x'| 

T = T + W, W{y) = j- f V ' T( *' } dV . 

4vr L, x-x' 
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The fields L, R, U, f , and T are by construction divergence-free. The vector potential tp of u satisfies 

VxV = ii, V-V> = 0. (16) 
For various purposes it is reasonable to consider dBi/dxj as constant. Then we have 

SB, aSjfx) r . , _ 1 f l.j(x') 3 



"«= 2 *j-£r- ^=sI^ dV - (17) 

We note that V • b = and V 2 b = — b. Analogous relations apply to V if dUi/dxj is constant. All 
quantities that occur in (|15p to Q17p may depend also on time, but for the sake of simplicity we have 
dropped the argument t everywhere. 



3 The traditional approach (approach (i)) to the mean electromotive force S 

In view of a later comparison with r-approaches we repeat here basic ideas of the traditional approach 
to the mean electromotive force £ and summarize some general results (see also, e.g., Krause and Radler 
1980 or Radler 2000). We recall that many applicat i ons h ave been made in the astrophysical context (e.g., 
Riidiger and Kitchatinov ( 19931 ) . iKitchatinov et al. ( 19941 )). For the sake of simplicity we exclude any mean 



motion of the fluid, that is, U = 0. If u is given, b is determined by equation (Jl]) and proper initial and 
boundary conditions. This equation is inhomogeneous in b due to the term with B. Any solution can be 
considered as a superposition of a solution of the homogeneous equation, which is independent of B, and 
a solution of the full equation which is linear and homogeneous in B. If there is a non-decaying solution 
of the homogeneous equation the mean electromotive force £ may have a non-decaying part, say 
independent of B (see, e.g., Radler 1976, 2000). If V x £^ = such solutions correspond to small-scale 
dynamos. Until further notice we assume however that b decays to zero as soon as B is equal to zero. In this 
case, discussed in section [37TT we speak of "purely hydrodynamic background turbulence". Later, in section 
3.21 we will admit a non-decaying b even if B is equal to zero and then speak of "magnetohydro dynamic 
background turbulence". We use the word "turbulence" here in a wide sense. If not specified otherwise 
u and b are considered as fields with arbitrary space and time dependencies but, of course, with zero 
averages. 



3.1 Purely hydrodynamic background turbulence 

3.1.1 Second order correlation approximation. We start here with the induction equation Q for the 
magnetic fluctuations b. As already mentioned we exclude any mean motion, that is, put U = 0. We use 
the second-order correlation approximation, here defined by neglecting the term G in equation @. Then 
this equation takes the simple form 

d t b - r]V 2 h = V x (u x B) . (18) 

Until further notice a non-zero divergence of u is admitted. We assume that equation (|18p applies in all 
infinite space. Then its general solution can be written as 

6 fc (x,t)= f G^(x-x',t-t )6fc(x',to)dV 

J oo 

+e k i m e m n P [ f G^\x-x\t-t')^(u n (x',t')B p (x.',t'))d 3 x'dt', (19) 
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where G^ is a Green's function, 

Gfii) (x, t) = G^ (x, t) = (4tt^)~ 3 / 2 exp(-x 2 /4r]t) , (20) 

and t = to means an initial time. We note that here, as a consequence of the neglect of G, solutions b 
which do not decay for B = are automatically excluded. 
For the calculation of £ we start from 

£i(x,t) = €ij k (uj(x,t)b k (x,t)} , (21) 

consider t — to as sufficiently large so that (uj(x, t) 6fe(x / , to)) = and let, for simplicity only, finally 
to — ► — oo. In this way we find 

£ i (x,t) = (e ijn 6 lp -e ijp 5 ln ) [ [°° r ) | Q jn ( x , f, -r) £ p (x -£,t-r) d 3 g dr (22) 

with the second-rank velocity correlation tensor Qij defined by 

Qi^t-^T) = {ui(x,t) Uj (x + £,t + T)). (23) 

We further assume sufficiently weak variations of B in x and t, and put under the integral in the sense of 
a truncated Taylor expansion 

5 p (x - €, t - r) = B p (x, t) - , (24) 

ignoring all other terms, that is, those with rdB p /dt, £ q !; r d 2 B p /dx q dx r , ... .In this way we find relation 
© with £^ = 0, that is £i = ai p B p + bi pq dB p /dx q , and 

dip = ~ [eijnOip - eijpdin) / / — — Qjn{x, t; £, -r) d 4 dr 

JooJo 5 

= (e ijn 6 lp - e ijp 5 ln ) [ [°° G^ (£, r) ^n(x^,-r) ^ ^ 

Joo JO 

i / r a \ f f°° 9G^r) ^ q 3 

JooJo <><, ? 



3.1.2 Range of validity. In order to formulate conditions under which the above results apply we 
introduce the magnetic Reynolds number Rm and the Strouhal number St, 

Rm = u c X c /r], St = u c T c /X c . (26) 

Here u c means a characteristic value of the fluctuating velocity u, and A c and r c are a characteristic length 
and a characteristic time of the variations of u. If u describes a turbulence A c and r c can be interpreted 
as correlation length and correlation time. We further define the dimensionless parameter q by 

q = X 2 J V r c . (27) 

It gives the ratio of the magnetic diffusion time \ 2 /r] for a fluid element whose size is characterized by A c 
to the time r c . We have then Rm/ St = q. If, e.g., St is fixed, the limit Rm — > coincides with q —* 0, and 
Rm — > oo with g — > oo. Sometimes the limits q — > oo and g — > are denoted by "high-conductivity limit" 
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and "low-conductivity limit" , respectively. This notation, introduced with the idea of varying rj for fixed 
finite values of A c and r c , might however be misleading. If, e.g., A c and 77 are fixed, q — ► 00 and q — ► 
correspond to r c —* or r c — > 00, that is, very quick or very slow variation of u in time. When discussing 
the limits q —* 00 or (7 — ► in the following we always consider A c as finite and t c as non-zero. 

The second-order correlation approximation applies as soon as the neglect of G in (j4]) is justified. A 
sufficient condition for that is |b| <C |B|. If q < 1 this is equivalent to Rm -C 1, if q > 1 to St <C 1. We 
may sum up this by saying that a sufficient condition reads 

min (Rm, St) < 1 (28) 

(see also, e.g., Krause and Radler 1971, 1980 or Radler 2000). Incidentally, determinations of and b^ 
have been done within the framework of numerical simulations of magnetoconvection or dynamos without 
using the second-order correlation approximation. It turned out that, e.g., for q < 1 there is a rather 
good agreement of the results for these quantities with those obtained in t he second-order corr e lation 
appr oximation not only for Rm <C 1 but for values of Rm up to about 10 ( Schrinner et al. 20051 . 2006 . 
20071 1. 



3.1.3 Limiting cases concerning q. Let us specify the above results for cii p and bi pq to the limiting 
cases with respect to q. As shown in appendix lAl in the limit q — ► 00 the relations (f25l) for ai p and bi pq 
turn into 

dip = (eijndip - e ijp 5i n ) [ (uj(x, t) dlln ^ t — L) dr (29) 

Jo oxi 

( f°° du n Cx,t — t) f°° \ 

= e ijnyj (uj(x,t) — — ^ }dT-5 np J (uj(x,t)(V • u(x,i - r))) drj 

/>oo 

bipq = -£ij P / (uj(x,t)u q (x,t - t)) dr . 
■> 

A contribution to bi pq containing 5 pq , which is because of V • B = meaningless for £, has been omitted. 
In the limit g->0we obtain 



£ijn$lp €ijp$ln f , , ,s , , /. .\\ > ^ £ 

— ^ — J (n i (x,t)n n (x + ^,t))6-^- 



{Trii K JjM^t) dUni Q^' t) ) y - 6np Ij u ^ t){v • u(x + tm T ) (30) 

l f 

bip q = ^ J (^ijnCp ~ tijpCn) £,q( u ji*-i t) U n (x + ^, t)} ^ . 

We note that these results for the li mits q — > 00 and q — > can also be derived from equation (I18p with 
r/V 2 b or 5^b, respectively, cancelled (jKrause and Radlerlll97ll . 1980l ). Interestingly enough, for the limit 



q — * 00 then, without introducing (1241) or similar assumptions, no contributions to £ with higher than 
first-order spatial derivatives of B occur. This is in agreement with the fact that, when considering first 
the general result for £ with higher spatial derivatives (Radler 1968a, b) and proceeding then to the limit 
q — > 00 all the coefficients of the higher than first-order spatial derivatives vanish. This finding might give 
some justification for considering no contributions to £ with higher than first-order spatial derivatives 
in astrophysical applications, in which the limit q — > 00 is a good approximation. There is, however, no 
general justification to ignore contributions to £ with time derivatives of B. By contrast, in the case q — > 
contributions to £ with higher spatial derivatives of B are well possible but no contributions with time 
derivatives of B. 
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3.1.4 Homogeneous isotropic turbulence. It is instructive to consider the special case of a homogeneous 
isotropic turbulence, for which £ takes the form ([7]), that is £ = aB — /3 V x B. For a. and /3, which are 
independent of position, we find then from (|25p 



a = ~\J J" 9G V d f T) ( u ( x > *) x u(x + t - r)) • (g/fl d 3 C dr 

= -1 / /°° (£, r) (u(x, • (V x u(x + £,t-r))) d 3 £ dr (31) 

Joo JO 

P = ~lLr ^If^ ^ *) • u(x + £, t - r)> £ d 3 e dr . 

Due to the assumed isotropy of the turbulence we may replace (u(x, i) • u(x + £, t — r)) by 3(u^(x, i) u^(x + 
t — r)), where U£ = (u ■ Incidentally, since G^ and all averages (...) depend on £ via £ only we 

may replace the integrals . . . d 3 £ by Air J* °° . . . £ 2 d£. 
In the limit g->oowe have 



i r°° 

a = - - J (u(x, i).(Vx u(x, t - r))) dr 



1 Z" 00 

/? = 3 y (u(x,t)-u(x,t-r))dr. (32) 

This result is often written in the form 

« = -i(u-(Vxu))r( Q ), /3=i(u 2 )r (/3) - (33) 

Here and r^) are just defined by equating the right-hand sides of (f32l) to those of (j33j) . respectively. 
This is, of course, only possible under the reasonable assumption that (u(x, t) ■ (V X u(x, t — r))) at r = 
is unequal to zero. It seems plausible to assume that both •n Q! ' and r^ are of the order of the correlation 
time r c of the turbulent velocity field. There is, however, hardly a general reason for their exact equality. 

Remarkably, as explained in appendix [HI in the case of an incompressible fluid and statistically steady 
turbulence the right-hand sides of (f32l) can also be expressed by integrals containing (u(x, t) ■ (V X u(x + 
£,t))) or (u(x, t) ■ u fx + £, t)} , taken over all £. 

In the limit q — > we obtain 

a = vt, £ (<u(x ' () x u(x + «• '» ' £) fr = " vb, l (u(x ' " ' (v x u(x + «' 4) » f 

"-sssjL^-^+^t- (34) 

Again the integrals . . . d 3 £ may be replaced by 47r J* °° . . . £ 2 d£. 

Incidentally, the results (|34p take a simple form if we represent ubyu = Vxi/> — V</> with the vector 
potential tp and a scalar potential 0. They read then 

a = -l(u-V) = -^(^-(VxV)), /3=-^«</> 2 >-<</> 2 » (35) 
077 377 077 

(see, e.g., Krause and Radler 1971, 1980, but note that there is a sign error in Krause and Radler 1971). 
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3.1.5 Comments. In the above calculation of the electromotive force £ no other restrictions concerning 
the fluctuating fluid velocity u have been used than some "smallness" which ensures the applicability of 
the second-order correlation approximation. 

We note that in this approximation if, as so far done, a purely hydrodynamic turbulence is considered 
and any influence of the initial b is ignored, £ does not contain a part £^ . 

We may of course specify the velocity u or its correlation tensor Qij such that they correspond to a 
turbulence on a rotating body, that is, subject to a Coriolis force. The dependence of u and Qij on the 
angular velocity ft responsible for the Coriolis force has to be calculated on the basis of the momentum 
balance (I14j) . Likewise we may consider a turbulence under the influence of a given mean magnetic field, 
that is, under the action of the Lorentz force, and specify u, again on the basis of the momentum balance, 
to be a function of B and its derivatives. In that sense the above derivations may serve also as a starting 
point for investigations of the quenching of the mean-field induction effects by the mean magnetic field. 

As for higher than second-order correlation approximations we refer to section [5TT1 

3.2 Magnetohydrodynamic background turbulence 

3.2.1 Second order correlation approximation. Let us relax the assumption that b decays to zero if 
B vanishes, that is, admit a magnetohydrodynamic background turbulence. We assume the fluid to be 
incompressible. In addition to the induction equation Q we use now from the very beginning also the 
momentum balance in the form (j!4[) . For the sake of simplicity we restrict our attention to a non-rotating 
fluid, that is fi = 0, and exclude again any mean motion, U = 0. Then u and b are governed by 

d t u = vV 2 u H ((B • V)b + (b • V)B + VM) + f + t, V • u = 

HQ 

<9 t b = ?7V 2 b + (B • V)u- (u-V)B + G, V-b = (36) 

with M, f and T as defined in () 15f> . 
We expand u and b with respect to B, 

u = u<°> + u« + . . . , b = b< > + b« + . . . , (37) 

where u*- -* and b^ ) are independent of B, further u^ 1 ) and of first order in B, and . . . stands for terms 
of higher order in B (see also Blackman and Field 1999 or Field et al. 1999). Then we have 

£ = £(») +£(!) + ... , fW =£(°1)+£(1°), ... 
£ (°) = ( U (°) x b(°)> , £^ = (u(°) x bW) , £W = ( U W x b(°)> , .... (38) 

We split now the equations (|36p in the usual way into equations for u(°) and b^ , which do not contain 
and b^, and equations for and bW, which contain of course also u*- -* and h^ '. We consider 
those for and b^ as being solved, that is, u(°) and b(°) as given. When deriving the equations for 
and b^ 1 ) we assume that f is independent of B. We further use the second-order approximation, 
in this context understood as the neglect of G^ 1 ) and in these equations. Here is defined by 
GW = V x ((u(°) x bW)' + (uW x b(°))') and f W is a quantity derived in the sense of CE3D and CE} from 
the analogously defined . In this way we arrive at 

tU m _ ^vV 1 ) = — ((B • V)b(°) + (b(°> • V)B + VM' 1 ') , V • u^ 1 ) = 

HQ 

d t b^ - r]V 2 b^ = (B • V)u(°) - (u(°) • V)B , V • b^ 1 ) = . (39) 



is defined as M in (TJS]) but with b(°) instead of b. 
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We recall here the derivations of section 13.1,11 Starting from equation (|18j) we have there calculated £ j 
in the form (|22j) . With the second equation (139p we find on the same way 



4 01) (x, t) = (e ijn S lp - e ijp 6 ln ) £ dGiV) d f T) ^ Qfl (x, t; -r) B p (x - t - r) d 3 £ dr . (40) 



Q^j is in analogy to Qij defined by 



Qif (x, t; r) = (^ (0) (x, t) uf (x + £ t + r)} . (41) 



Because of V • u<°) = we have dQ^f/dxi - dQ^f/d^ = and dQ^f/d^ = 0. 
The first equation (|39|) leads to 



f?^fivf(x, t ;- £ ,- T )d^. (42) 



, r) is defined like G^(£ , r) in (I20j) . but with r/ replaced by za Further - is in analogy to Q 
defined by 



'J 



i?^(x,t;C,r) = (6W(x,t)6W(x + €,i + r)) (43) 



A^(x t-£ t)- — [ a <W;*',T) gg^x+j^+r) 



and ivj 0) by 



Because of V • b<°) = we have dR^/dxi - dB$ /d£i = and dllf) /d^ = 0. 

We restrict our attention now to the limit of small B. As in section [3. 1.1 1 we assume in addition sufficiently 
weak variations of B in x and t and introduce the truncated Taylor expansion (|24p into the integrands of 
(1401) . (|42f) and (I44p . Considering ([38]) and ignoring contributions to £ of higher than first order in B we 

find then again ©, that is £i = £^ + a^- -B-, + bij^dBj/dxk- Here we have = (l/ ) x b^), and this 
may in general well be unequal to zero. We put 

a- ■ - a {u) + a (b) b- u - b {u) + b (b) (45) 



with a£j and 6^ determined by u^°\ and afj and bf^ k by b(°\ Then we have 



oo JO 



3 

(0), 



Joe JO 



p Joe Jo °Z,P 
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« ~*" £ f^Pf - T )# £dT (46) 



^ipq 



J G<"> (f , t)Q(°» (x, t; { , -r) d 3 f dr 

/ / GM (£, r )(ey P - 2 e 07 — — Aj 1 )^ (x, t; £, -r) d 3 £ dr , 

7 oo 70 OQO^p 5 



W Joo JO 



A ? F(x,^,r) = --y oo df. 



(47) 



By the way, the above relations for af? and b^ q can also be concluded from (j25l) by replacing a^, 6j M and 
Q jn with a-p\ 6^ and Q^, respectively, and considering dQ^'/d^j = 0. 



3.2.2 Range of validity. In addition to the dimensionless parameters Rm, St and q introduced with 
(j26j) and (|27l) we define the hydrodynamic Reynolds number i?e and a parameter p by 

Re = u c \ c /v , p = \ 2 c /vt c . (48) 

We might call the cases p — > oo and p — > "low-viscosity limit" and "high-viscosity limit" . However, the 
remarks made on the notations "high-conductivity limit" and "low-conductivity limit" made under (I27j) 
apply analogously. As in the case of q, when discussing the limits p — > oo or p — > in the following we 
always consider A c as finite and r c as non-zero. 

Modifying the reasoning which lead us to the condition (128ft properly we find that a sufficient condition 
for the validity of our above results for magnetohydrodynamic background turbulence is given by 

min(Rm, St)<l, min(i?e, St) < 1 . (49) 



3.2.3 A limiting case. We note that and 6^ depend only on q and not on p, and af^ and bf^ k only 

(u) (u) 

on p and not on q. The relations for a- and 6^ in the limits q —* oo and g — > agree with those for Oy 

and ftjjfc given with (J29I) and (i30j) . specified to V • u = and with u is replaced by u^ ). For af^ and bf) k 
we find on the way described in appendix [X] in the limit p —* oo 

°- = " wio ^ (X ' } } 



/>oo 

6< 6 > 

tpq 



(5 0) (x,t)6f (x,t-r))dr (50) 

HQ Jo 
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A contribution to of^ q containing S pg has been omitted. In the limit p — > we have 

C = -4^ jjS^, + i m Jf-){bf (x,t)6<?>(x + «,t)>^ (51) 



3.2.4 Homogeneous isotropic turbulence. Let us now again consider the special case of homogeneous 
isotropic magnetohydrodynamic turbulence. Then £ takes again the form ([7]), that is £ = aB — /3 V x B. 
On the basis of (SSI) and (l4"6l) we find 



a = a 



M - a (b) , (3 = pW> (52) 



with 



= -i 11°° T) (u(0) (x ' ^ x u(0) (x + * - r)) • m) d ^ dT 

= _H / OO GW(e,r)(u( )(x,t).(Vxu(°)(x + ^,t-T)))d 3 edr 

Joo JO 

i /■ rsc^ (0) (0) )t _ 3(dT (53) 

Joo JO 

- / / OO GW(C,r)(b( )(x,t).(Vxb(°)(x + ^t-r)))d 3 edr 
W Joe Jo 

= I / /°°GW(e,T)(u( )(x,t) •u(°)(x + ^,t-r))d 3 edr. 

3 Joo JO 



a< 6 > 



3w. 



By the reason explained in the context of (|31|) the integrals f . . . d 3 £ in (|53j) may be replaced by 
47T f °° . . . £ 2 d£. Interesting l y eno ugh there is no "magnetic part" of (5. This agrees with a result by 
Vainshtein and Kichatinovl ( 19831 ) obtained within an early version of the T-approach applying to the 



limit q — > oo. 
In the case q,p — > oo we have 



1 f 00 

■- / (u (0) (x,t).(Vxu<°)(x,i-T)))dT 

3 Jo 

a w = / (b(°> (x, t) • (V X b^ (x, t - r)))dr (54) 

3/U£» Jo 



1 /'°° 

= - / ( u (°) (x, t) ■ u(°) (x, t - r))dr . 
3 Jo 



In the sense discussed in the context of (1331) we write this also in the form 



a H = _I/ U (0) . ( V x u(°))) r( Q «) , = -— (b(°) • (V x b(°>)> 
3 3/j.q 

P<M) = I( u (°) 2 )r (/3n) (55) 
3 
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with times and r^ u ) defined by equating the corresponding right-hand sides of (|54p and (|55j) . 

3.2.5 Comments. In the results for purely hydrodynamic background turbulence presented in section ETTl 
the velocity u may well depend on B. If we ignore this dependency u agrees with u(°) introduced here. 
In this case the results of section [3T1 coincide with those obtained here if specified to h>(°) = 0. We stress 
that in the case of a purely hydrodynamic background turbulence, even if we assume that u is subject 
to an influence of the Lorentz force and so depends on B, e.g., the relations (f25l) . (|29j) and (i30j) for ai p 
remain their validity and have not to be modified by adding terms which explicitly depend on b. Only if 
we replace there u by u^ ^ an additional term containing b might occur. 

The relation for a defined by (f52l) and (f55|) is usually taken to be supported by the analysis of 
Pouquet et oil (jl97fih . who employ a much more sophisticated closure than the r-approaches, that is, 



the "eddy damped quasi-normal Markovian approximation" . We note however that the justification given 
there (subsection "Phenomenology of the M-helicity effect" of their section 3) is based on the assumption 
of an (at least initially) purely magnetic background turbulence, see their equation (3.19). It has to be 
considered as a consequence of a given electromotive force rather than a kinetic driving force. Under these 
conditions, too, it is not surprising to see a term depending on the current helicity (b • (V x b)) in the 
expression for a, just as in our above derivation leading to the results (|52p and (|55p . 



4 The t approaches (approaches (ii)) to the mean electromotive force £ 

As mentioned above, several analytical approaches to the mean electromotive force £ going beyond the 
second-order correlation approximation have been propose d which mak e use of some modification of the 
T-approximation of the hydrodynamic turbulence theory (jOrszag 1970I ). In contrast to approach (i), in 



these approaches, labelled by (ii) in the following, both the induction equation @ and the momentum 
balance in the form (|14p with the forcing term f and the non-linear terms G and T are used from the very 
beginning. In this way also V • u = is introduced. Equations ([!]) and (|14p are however not really solved. 
Instead, relations for £ are derived which contain, of course, contributions resulting from the mentioned 
quantities, that is, f , G and T. For these contributions then some kind of T-approximation is introduced. 
In this section we give a critical review of the approaches of type (ii). In particular we check to which 
extent their results satisfy the elementary requirement that they agree with results of approach (i) in the 
range of its validity. 

For the sake of simplicity we exclude again any mean motion of the fluid, that is, U = 0. We further 
restrict our attention to the case of a non-rotating fluid, that is, $7 = 0, and give only short comments on 
the case with ft ^ 0. 



4.1 The simple r approach (approach (iia)) 

Let us first describe a simple approach to the mean electromotive force £, which we call approach (iia) 
in the following. It defines a frame which allows us to e x plain and to discuss the approaches used 
by Vainshtein and Kichatinov ( 19831 ). Blackman and Field! |2002h and iBrandenburg and Subramanianl 
(|20n5bh . 



4.1.1 Outline of approach (iia). We calculate first the time derivative of £ on the basis of 

d t £ = (d t u x b + u x d t h) . (56) 

Using equations further (fT4"|) with R = U = 0, and considering dB p /dx q as constant so that (fT7|) 
applies, we obtain 



8 t £ = I + Y + Z , (57) 
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where 

du k 1 db k 1 d^b dB 

h = £ijk{(uj^—} ( b j^—)) B p ~ £ijk((ujU q )6 kp ((b k b q }5 jp + 2 (- — b k )))—-Z 

OXp HQ OX p fJ,Q OXjOXp OX q 

y% = £ijk(v(Auj) b k + r/UjAbk) + e ijk (fj b k ) (58) 
Zi = £ijk{Tj b k + Uj G k ) . 

The term Y depends on v and rj. Even if these quantities vanish it remains non-zero as long as €ij k {fj b k ) 
does so. The term Z considers the non-linear terms G and T that occur in @ and f|14)) . It vanishes if 
these terms do so. 

Let us now split £, I, Y and Z into parts £^°\ Y^ and Z", which are independent of B, and 
remaining parts £^ , 1^ B \ and Z^ B l Clearly = 0, that is = I. We do not deal with f (°) 

and remark only that it is equal to zero for several simple cases. For £^ we have 

d t £^=I + Y^ + Z^. (59) 

We introduce now the assumptions 

Y^ = -£^/t y , Z^ = -£^/t z (60) 

with two times Ty and t#, which are considered as non-negative. The second of these equations introduces 
some kind of r-approximation, which reduces terms of third order in u or b to such of second order. 
When accepting (|60p we may rewrite (|59p into 

dtEW + (l/r«.)£ (B) = I , l/n = 1/ty + l/r z . (61) 

This is equivalent to 

£i B )(t) = £ {B) (t ) exp(-(t - t )/n) + f ° I(t - t') exp(-tVn) dt' , (62) 

J o 

where to means some initial time. We let to — > — oo and assume that I does not markedly vary in time 
intervals with a length of some r*. In this way we find a result for £ which we can write in the form of (jGj) , 
that is £i = £^ + ai p B p + bi pq dB p /dx q , with 

(I duk \ 1 lu dhk \\ 
a ip = e ijk (( Uj — )--{b^ »r, 

1 d 2 b 

h pq = -eij k {(uju g }5 kp ((b k bq)5 jp + 2 ( 9 b k ))) r* . (63) 

f-LQ C/X j C/Xp 

Let us again consider the special case of a homogeneous isotropic turbulence. Then £^ has to be equal 
to zero and we have again Oj p = a5i p and bi pq = (3 ei pq . Consequently £ takes the form (J7J), that is 
£ = aB-/?VxB. With ([63]) and considering V 2 b = -b we find 

a = ~((u-(Vxu))-i (b'(Vxb)»r„ /3=I(n 2 ) n . (64) 
3 hq 3 

There is a kinetic and a magnetic contribution to a but again only a kinetic one to (5. 

It is of some interest to consider also the case of a homogeneous turbulence deviating from isotropy 
only due to the presence of a homogeneous mean magnetic field B. Then we have ai p = a\ Si p + c^eejeBp, 
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where e B = B/|B|. This leads to £ = [a.\ + a2)B and a.\ + «2 = ai p e B je Bp . We may write again £ = aB. 
However a has then no longer the form given by (1641) but 

a = ((u x (e B • V) u) • e B - — (b x (e B • V) b) • e B ) r* . (65) 

It is in any case the dependence of u and b on B which leads to a dependence of a on |B|. In addition a 
dependence of r* on |B| cannot be ruled out. 

In the limit of small B the averages (...) in (|63|) and (|64|) turn into (. . .)(°), what means that u and b are 
replaced by the corresponding quantities u.(°) and b^ ) for B — > 0. In view of this limit we further remark 
that the average of the three values of a given by ([65]) with e B = (1, 0, 0), e B = (0, 1, 0) and e B = (0, 0, 1) 
agrees just with a given by (|64p . again with u and b replaced by u(°) and b(°). 

4.1.2 Comments. In the above derivation of the result (|63p the assumptions (|60p play a crucial role. 
If we simplify our reasoning by considering (|6ip from the very beginning for the steady case only, a very 
strange aspect of the way to a relation of type ©, £j = . . . + ai p Bp + bi pq dBp/dx q , becomes obvious. There 
is no original relation of this type, which then would be improved by the assumptions (|60p . Instead, this 
relation is just established by the assumptions (|60p . £ on its left-hand side originates only from the £^ 
in (f50|). 

Nevertheless, at the first glance the assumptions ()60|) look plausible at least for simple cases. For example, 
for an isotropic background turbulence and a homogeneous steady mean magnetic field B the quantities 
Y( B ), Z( b ) and £^ must by symmetry reasons be proportional to B. This leads immediately to relations 
like (|60p . The minus signs in these relations together with positive values ry and tz ensures that £ remains 
finite. 

The assumptions (|60p become however questionable in more complex cases. Consider, for example, again 
an isotropic background turbulence but admit a inhomogeneous steady field B. Clearly Y^ B \ and 
£( B ) are determined by the vectors B and V x B and must have the structures c\ B + C2 V x B with two 
scalar coefficients c u , v = 1, 2. However, assumptions (I60|) could only be justified if the ratios jci^ 1 and 

(ffl l<ffl of the coefficients c„ , c„ and c„ for Y( B \ Z^- 8 ) and £^ B \ respectively, were independent of 
v. We see no reason for that. Some way-out could consist in splitting ~Y^ B \ and £^ in terms which 
correspond either to c\ B or to C2 V x B and to formulate ansatzes like (|60p for each term separately. 

Even if we accept the assumptions ([60]) the question remains about realistic values of ry and tz, that 
is, which value of r* should be inserted in the results. 

Let us point out a conflict of the results of approach (iia) with those of approach (i). Consider the special 
case in which the nonlinear terms G and T tend to zero, that is, Z — > 0, or tz — > oo. We have to require 
that then results of approach (iia) like (|63p or (|64p turn into the corresponding results of approach (i). 
There is indeed some similarity of the results (|63p and (164ft with the results (129ft and (|33p . or (I50p and (1551) 
of approach (i) in the limit q,p — > oo . However, (]63p and (164ft do not reflect the fact indicated by (I29p 
and (I50p . or (|32H and (|54f) . that ajj and bijk, or a and /3, in this limit depend on correlations between the 
components of u at two different times. According to the above reasoning the results (|63p and (164ft should 
at the same time correspond to the results of approach (i) for the limit q,p — ► 0, that is, (f30l) and (|34l) . or 
(|5Tj) . This applies obviously at best in a very crude sense. In particular, the dependence of aij and bijk, or 
a and /?, on correlations of components of u in different points in space is not reproduced. 

There is another remarkable difference of the results of approaches (iia) and (i). Consider the case of a 
purely hydrodynamic background turbulence. In approach (i) the quantities and bijk as given by (f25l) . 
()29p or (I30p depend only on the correlation properties of u but, even if u is influenced by the magnetic 
field, not explicitly on those of b; see also section f3.2.51 In the relation (I63p of approach (iia), however, 
in addition to terms with u also such with b occur, which need not to be negligible beyond the limit of 
small B. In the relations for a given with and ([65]) terms with b occur, too. Whereas (f64"|) applies 
to isotropic turbulence and thus seems reasonable for small B only, (I65p can be considered for larger B, 
where this term needs no longer to be small. 
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Admitting again magnetohydro dynamic background turbulence we further note that the times T^ au ^ 
and 

7>&) in 

relation (|55j) for a, which has been derived in approach (i), need not to coincide. In the 
corresponding result (|64h of approach (iia), however, only one comparable time, r*, occurs. 

Since there is hardly any reason to doubt in the results of approach (i) we conclude that those of approach 
(iia) are, at least in the range of the validity of approach (i), not completely correct. This, of course, calls 
also results of approach (iia) for more general cases into question. 

We further remark that there is no straightforward extension of approach (iia) as presented above to 
the case of turbulence in a rotating system. In this case the equation that occurs instead of (I57p contains 
in addition to £, which is defined by a product of u and b with both taken at the same point, also a term 
with a product of u and b taken at different points. Then £ can no longer be determined on the simple 
way used above. It is to be expected that aij and 6^ contain correlations of components of u at different 
points in space. 

Finally, even an extension in that sense can hardly deliver a correct result for a contribution to £ like 
-<5i(« • V)B in ©, because of (ft • V)B = — ft x (V x B) + V(fi • B) often discussed as "ft x J-effect". 
In approach (i) this contribution results, at least in the case of homogeneous turbulence, from a part of 
the correlation tensor Qij of u which is odd in the difference r of the times at which the components of 
u are taken, see appendix This part, however, does not occur at all in approach (iia), which considers 
only correlations of the components of u at the same time. 

The relation (1641) for a is often used as a starting point for investiga tions on a-quenching, e.g., in 



Blackman and Field! (|2002h and in IBrandenburg and Subramanianl (|2005al lbl) . We point out that this re- 



lation applies only for isotropic turbulence, that is, in the limit B — > 0. Investigations on a-quenching 
carried out within the framework of approach (iia) should rather start from a relation like ()65|) . On this 
level, of course, instead of the simple relations for (b • (V x b)) used in the mentioned investigations, 
which have been concluded from the magnetic helicity balance, corresponding relations for quantities like 
(b x (ee • V) b) • ee are needed. Moreover, as we know from the above comparison with approach (i), in 
the case of a purely hydrodynamic background turbulence the existence of a magnetic contribution to a 
and so its role in a-quenching is questionable. 



4.2 A spectral r approach (approach (iib)) 

Let us consi der the approach used in the papers by Radler et al. ( 20031 ) a nd by Rogachevskii and Kleeorin 
(|200l l2004!l. It has been extensively repeated in the survey article by IBrandenburg and Subramanianl 
(20053), who name it "minimal r-approximation (MTA)". As mentioned above we exclude here again any 
mean motion and restrict ourselves to a non-rotating fluid, that is U = ft = 0. 



4.2.1 Outline of approach (iib). In view of the determination of £ the attention is first focused on the 
cross-correlation tensor Xjk defined by 

Xii(x,0 = <«i(x + €/2)6i(x-C/2)). (66) 

Here and in what follows all quantities are considered at the same time. For the sake of simplicity we drop 
the arguments t without ignoring the dependence on t. 
In addition to Xij other correlation tensors are of interest, which we first define in the general form 

%(v, w; x, = ( Wi (x + £/2) Wj (x - £/2)> . (67) 

We use a Fourier transformation with respect to the space coordinates only, F(x) = J F(k) exp(ik-x) d 3 k. 
As explained in appendix [D] the Fourier transform <3?jj of <&ij with respect to £ can be represented in the 
form 

$ 4i (v,w;x,k) = f(vi(k + K/2)wj(-k + K/2)) exp(iK • x) d 3 if . (68) 
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For the calculation of £ it is sufficient to know the antisymmetric part of the tensor Xiji that is, the 
vector x defined by Xi = e ijk Xjk- We introduce its Fourier transform with respect to 

Xi(x, k) = e ijk Xjjfc(x, k) . (69) 

Then we have 

£{x) = J x(x,k)d 3 fc. (70) 
Let us consider the time derivative of Xij, 

Qtxy(x,k) = J (^(k + K^jSj-C-k + K^+UiCk + K^J^S^-k + K^)) exp(iK • x) d 3 K . (71) 

We consider £ and the quantities necessary for its determination at a given point and choose the coordinate 
system such that there x = 0. In that sense we assume for the calculation of d%U{ and dtbj simply 

Bi(x,t) = B t + B ijXj (72) 

with constant Bi and Bu. Then we obtain with the equations (j3J), (|14|) and (I15h . reduced to CI = U = 0, 

1 / - - kki- <%fkK 

dtUjQt) = -ukHjQt) + —[WfabjQt) + Bjibi(k) - 2B lrn ^-b m {k) - B lm k-£-t) 

+/ i (k)+T J -(k) 

d t b k (k) = - V k 2 b k (k) + LB,fc { «fc(k) - S««,(k) - Amfci^ + G fc (k) . (73) 



"in 



We insert this into the integral (|71[) and split it then into a sum of integrals of the type 

J F(k, K) (uj(k + K/2) Wjt(k — K/2)) exp(iK • x) d 3 ^ . (74) 

We may expand F in a power series with respect to K and rewrite each factor if; into a derivative 
Vi = d/dxi of the integral. Since we intend to ignore all contributions to £ with higher than first-order 
spatial derivatives of mean quantities we do the same with respect to Xij- I n this context, of course, Bn 
has to be considered as derivative of the mean quantity B{. Using then (f69j) we obtain 



with 



dtX = -{y + r])k 2 x + \{y - r?)(k • V) X + I + Y + Z . (75) 



h = tijk ( ~ iBikiu jk (x,k) + ^BiViii jk (x,k) - B kl Uji(x,k) 
-\Bi m ki dUjl f^' — - \Bi m kix m u jk (x, k) 



1 / 1- - kki- 

+— [iBikib jk (x, k) + ^BiVi b jk (x, k) + B jZ 6 Zfc (x, k) - 2Bi m -j^-b mk (x, k) 



~Bi m ki ^q*' k + iBi m kix m b jk (x, k) ) j M>) 
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% = e ijk $ jk (i,b;yL,k) 

Zi = e ijk (* jfc (f , b; x, k) + * ifc (u, G; x, k)) . 

The Ujk and bj k are the Fourier transforms of Uj k and bj k defined in the sense of ([671) by 

%Jfc(x,0 = $,-fc(u,u;x,£) , &jfc(x,£) = $j fe (b,b;x,|). (77) 

In contrast to Y used in approach (iia), Y no longer contains the diffusion terms, that is neither v nor r], 
but is only a consequence of the forcing term in the momentum balance. Like Z in approach (iia) also Z 
describes the influence of triple correlations of u and b on £. 

We split now £ and likewise y, i, Y and Z into parts £^\ X^°\ Y^ -* and Z^ ) independent of B 
and remaining parts £^ B \ x^ B \ 1^ B \ and 7S B \ We have, of course, = 0, that is, 1^ = 1. Again 
we do not deal with £^ . In view of£( s ) we conclude from (1751) that 



d t ± {B) = -{y + v)k 2 X (B) + K v ~v)(k- V) X {B) + I + Y (B) + Z< B > . (78) 
In analogy to ([601) we introduce here the assumptions 

Y( fi )(x,k) = - X {B) (x,V/f Y (k), Z( fi )(x,k) = -x (B) (x,k)/f,(fc) (79) 



with positive functions fy and fz of k. In Radler et al. ( 20031 ) it has not been distinguished between the 



two terms Y^) and iS- B \ although they are of quite different origin. For their sum an ansatz like that for 
Y( fi ) or Z( fi ) is used. This is, however, of minor interest for the following discussions. 
Starting from ([78D and using (179p we obtain 



0t X (B) + (Vn )X {B) = i (f - r/) (k • V) X {B) + 1 , 1/f* = + r?) /c 2 + 1/fy + l/f z . (80) 

We recall that all contributions to £ with higher than first-order derivatives of mean quantities should be 
ignored. This applies, of course, to x^ and i, too. In that sense we conclude from (|80l) that 

d t (k • V) x (B) + (Vr*)(k • V) x {B) = (k • V) i . (81) 

In view of £ at x = we consider also the equations (|8Q|) and (j8ip at x = 0. They are then differential 
equations for x and (k • V)x with respect to t. Solving them in the same way as we did it with (j61h . 
assuming here, too, to ~~ * ~~ °° and sufficiently small variations of I in time, we find 

X (B) = (i + i(^-r ? )(k-V)in)n, (82) 



and, together with ([76]). 



1 11- 



+(u - r)) kk p Vi (u jk - — b jk ) %(k)) f*(k) B p 
- (u jq 5 kp + ^k p -^u jk - (b qk 5 jp -^k p -^rhk ~ 2 b qk ) (83) 
-(u-rj) kpk q (u jk — bj k ) f*(A;)) T*(k)B p 



"pq J j 



k 
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again with all quantities taken at x = 0. 
With this result we can calculate £ according to (I70p . Admitting again arbitrary x and writing the result 

in the form ([6]), that is, Si = S^ + a%jBj + bijkdBj/dxk, we have 

aip = e ijk J Qv p u jk {x, k) - ik p Uj k (x, k) + b jk (pc, k) + ik p b jk (x, k) 

+{v-r{) kik p Vi ^n jfc (x,k) - ^- b jk (x, k)) f*(£;)^ f*(fc)d 

frjpg = ~£ijk J ^jq( x > k )$fcp + — k p -^j—Uj k (x., k) 

-(i/ - 77) k p k q (uj k (x,k) - — 6 jfc (x,k)] f*(/s)) f*(/c)d 3 fc. 
v fig / / 

These relations are by the convolution theorem equivalent to 

€i i k ( f (l ( , tin + A^( X ~ ^/ 2 '*)\ 



dx p 



+1 _ B bj (, + (/2,t) b ]r 

\1Q OXp I 



— (bj( X + ^/2,t)b k ( X -^/2,t))) r-(i)d-\ 

HQ 



— ((6,(x + £/2,t) & fc (x-£/2,t))J ip 



<M X + €/2, t) M x - Z/2, t)) 



2 

2 ^-(A-^b^x + £/2, t)6 fc (x - £/2, t)» ) ) r*(£) d 3 £ 

- — (&,(x + 4/2, t) 6 fc (x - 4/2, t)>) r, 2 (0 d 3 £ ] . 

For clarity we have inserted again the arguments t. The quantities r* and r* are given by 

n(k) exp(ik • 4) d 3 Ac = — J t*(k) sm(ktl)kdk 



(84) 



b m = -J^}( / ((%( x + ^/2,t)^( x -4/2,t))4 P 

(n j (x + 4/2,t)n fc (x-4/2,t)) (85) 
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rH0 = J T* 2 (A0e X p(ik.Od 3 A; = ^3 £ r*(|£ - r*(£') d 3 £' . (86) 

Note that r* has the dimension time/length 3 , and r* the dimension time/length 3 ^ 2 . As for A^ 1 we refer 
to the definition (gZD. 

Let us again consider the special case of homogeneous isotropic turbulence. Then we have again £^ = 
and dip = ci 5i p and bi pq = (5 6i pq so that £ takes the form ((ZJ), that is £ = a B — (3V x B. Using (i85j) we 
find 

« = -37^3 / (< u ( x + ^/ 2 ' *) • ( V x u ( x " £/ 2 ' '))> 



— <b(x + i/2, t) ■ (V x b(x - i/2, t))j) r*(£) d 3 £ 



/9 = 3^3 £< u ( x + *A *) • u ( x " €/2, 0) MO d 3 £ . (87) 

We have here replaced (b(x — £/2,i) ■ (V x b(x + £/2,i))), what originally occurs in the derivation from 
(1851) . by (b(x + £/2, i) • (V x b(x — £/2, i))). Due to the assumed homogeneity of the turbulence the averages 
(. . .) in (|87p do not depend on x , and we may further replace, e.g., (u(x + £/2, t) ■ (V x u(x — £/2,t))) with 
(u(x, t) ■ (V x u(x + i))). Again, by the reason explained in the context of (j3Tj) the integrals . . . d 3 £ 
may be replaced with 4-7T J °° . . . £ 2 d£. The — 77) terms in (1851) are without influence on a and /?. 

In Radler et al. (j2003h the possibility of a non-zero has been ignored from the very beginning. In 



addition only the limit of small B is considered. So for a comparison Uij and bij in (|83p and (|84p have to 
be interpreted as and , the limits of these quantities for B — > 0, and likewise u and b in ([85]) and 
(j87[) as u(°) and , understood in the same sense. 

4.2.2 Comments. Similar to the situation with approach (iia) and the assumptions (|60p on Y^' and 

the results of the above derivation depend crucially on the assumptions (179p on Y^ and 
describing the influences of the forcing and of the triple correlations of u and b. So we have here essentially 
to repeat the comment made in 14.1.21 in view of the assumptions (|60p used in approach (iia). Here, too, 
we see no rigorous way to justify assumption (|79p . Moreover, the choice of fy and fz remains open. 

A very strange aspect of the above procedure is that it needs as input no other correlation tensors than 
those with components of u or b taken at the same time. 

Let us again consider the case in which the nonlinear terms G and T tend to zero. This implies 0, 
or tz — > 00. In this case the results of approach (iib) have to coincide with those of approach (i). 

As can be seen, e.g., from (|80p the results of approach (iib) depend in general on both v and r/. This 
is clearly in conflict with the fact that those of approach (i) for purely hydrodynamic turbulence do not 
depend on v. 

In the limit q — > 00 approach (i) delivered us the quantities a%j and in the form (|29p . These are far 
from coinciding with the results (|85p obtained in approach (iib). A striking discrepancy is that in approach 
(i), apart from the independence of the results on v, only correlations of the components of u at the same 
point in space but at different times occur, in approach (iib) only correlations at different points in space 
but at the same time. At the first glance this discrepancy seems to disappear in the case of homogeneous 
isotropic statistical steady turbulence, in which the results of approach (i) for a and (5 given by (i32|) can 
be rewritten in the forms (|B9p and (jB12[) . which correspond to ([87]) . While however the r c (IB9P and (]B12[) 
may differ from each other the r* in (j32[) have to coincide. 

In the limit q — > approach (i) delivers us the relations (|30l) for aj,- and bijk- They are in so far similar 
to (I85p as in both cases only correlations of the components of u at the same time occur. Nevertheless a 



coincidence of the two kinds of results would require very special properties of the factor r» in (i87|). We see 



no reasons for such properties. When putting on trial Y( s ) = ( Brandenburg and Subramanian ( 2005al ) 
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consider ~Y~( B ^ as negligible) we arrive at expressions for a and f3 which differ from (|34p only in so far as 
rj + v occurs instead of rj. A view at the first of the equations (fT3j) might suggest to put = vk 2 x . 
Then we obtain indeed a and (3 in full agreement with (|34p . Nevertheless = vk 2 x^ removes v only 
from one of two places at which it occurs in (|80p , and even if it does no longer occur in a and (5 it will do 
so in other mean-field coefficients. We believe that an essential reason for shortcomings of approach (iib) 
is a treatment of the term Y^ B \ which is in general not justified. 

In the discussion of the results of approach (iia) we have pointed out a difference to those of approach (i) , 
which becomes visible in the special case of purely hydrodynamic background turbulence. This difference 
occurs here, too. In this case in approach (i) the quantities aij and 6^, or a and /?, can always be 
represented by the correlation properties of u. In approach (iib) however, beyond the limit of small B, the 
expressions for these quantities except (3 contain both u and b. 

We see that approach (iib) to £ delivers results which even under comparable conditions deviate from 
the results of approach (i), which we consider as correct. Therefore also beyond the range of validity 
of approach (i) and beyond the case of non-rotating turbulence, all results derived in the framework of 
approach (iib) have to be considered with caution. Their details are in general not correct. 

If we admit a Coriolis force, that is ft ^ 0, or a mean motion, there is a coupling between the 

symmetric and antisymmetic parts of and the calculation of x is more complex. By the reason already 
discussed in section H. 1.21 it is then hard to imagine that an extension of approach (iib) to the case ft ^ 
can describe, e.g., the 11 ft x J-effect" correctly. Here a view on the "ft x J-dynamo" is of some interest. 
Its existence depends on a non-zero value of a specific coefficient. In Radler et al\ ( 20031 ) . equations (130) 
and (132), it is given by 5q — kq, and for hydro dynamic background turbule nce it is non-zero as long as f* 
depends on k. According to the calculations hv kadler and Stenanovl » done again for hydrodynamic 
background turbulence and in the second-order correlation approximation, however, this value vanishes in 
both the limits q — ► and q — ► oo independent of any assumption on the correlation time; see section VD 
a nd appendix C of that pape r. 

iRogachevskii and KleeorinI jjjoal ). considering a mean motion and proceeding in the spirit of approach 
(iib), found a contribution to £ proportional to W x (V x B), where W stand s for V x U. A contribution 
which can be interpreted in that sense occurs already in a paper bv lUrpinl (|l999T l. though as a result 
of a not easily comprehensible calculation. The existence of thi s "W x J-effect" has been confirmed in 
the framework of the second-order correlati on approximation bvlRadler and Stepanovl (|2006h . The simple 
model of an "W x J-dynamo" proposed by Rogachevskii and Kleeorin ( 20031 ) (see also Rogachevskii and 
Kleeorin 2004) works with the coefficients determining this and the accompanying effects as calculated in 
the spectral T-approximation. It turned out however that it fails to work with the coefficients resulting 
from the second-order correlati on approximation for purely hydrodynamic background t urbulence, and 
this applies independently of q ( Riidiger and Kitchatinov 20061 . Radler and Stepanov 20061 ). 



4.3 The t approaches in the light of numerical simulations 

Quite naturally the question arises whether the correctness or t he limits of validity of the r-approaches 
can be checked by means of direct numerical simulations. Indeed, Brandenburg and Subramanian ( 2005bl ) 



report on numerical experiments with forced helical magnetohydrodynamic turbulence in a rectangular 
box aimed at studying the a-effect and a-quenching in the framework of approach (iia). In the underlying 
theory they use the ansatz (|60f) for (in their notation T), standing for triple correlations, but ignore 
~y( B \ which comprises the effects of viscous and ohmic dissipation and the correlation of the magnetic 
fluctuations with the external force. 

Brandenburg and Subramanian examined the relation between and a dynamo-generated mean 

field B on the basis of a simulation with a forcing yielding an isotropic background turbulence. Under the 
assumption of the proportionality of £ and B, justified by a weak spatial variation of B, the result turned 
out to be consistent with the ansatz (160p for 7S- B \ The simulation confirmed not only the proportionality of 
Z( B ) and £, which was due to the isotropy of the turbulence to be expected already by symmetry reasons, 
but delivered also a non-negative tz- The ansatz (160p for Y^, which was not considered in this work, 
is plausible by symmetry reasons, too. So under the considered circumstances the neglect of a non-zero 
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Y( B ) cannot have any other consequence in the result for a than a deviation of the time r* from rz- 

In a further step Brandenburg and Subramanian extracted the time r* (in their notation r) from the 
simulations by making use of the relation (164p for a. Of course, using (|64h instead of (I65h has generally to 
be questioned as only the latter is the relevant one beyond the limit B — » 0. The use of (164p can only be 
justified under the additional assumption that — 3 (u • (V X u)) and -| (b ■ (V X b)) are proper proxies of 
(u x (ee • V) u) -ee and (b x (es ■ V) b) -ee- In some cases, from the simulations negative values of r* emerged 
(surprisingly enough just for small B), what is in contradiction with the idea of a mean electromotive force 
decaying in the absence of a mean magnetic field. As a way-out an ad-hoc modification of (164p was 
introduced in form of additional quenching functions, gx and g^, in front of both the kinetic and the 
magnetic contributions to a. This measure cannot be brought into accordance with the derivation of (|64h . 
which allows at best the same quenching function for both contributions, that is, a dependence of on |B|. 
Interestingly enough, in the relations (|52p and (|55p for a, derived in approach (i) for magnetohydrodynamic 
background turbulence in the limit q,p — > 00, the two times r( an ) and r( afe ) occur, which could be well 
different from each other. As far as the assumption of magnetohydrodynamic background turbulence, that 
is, of a small-scale dynamo, is justified, the numerical results support our doubt in the correctness of 
approach (iia) based on the conflict of the latter with approach (i). It is however not clear to us whether 
all runs correspond to this assumption. In the case of a purely hydrodynamic background turbulence there 
is still another doubt in relation (|64p for a, which results also from the comparison with approach (i) and 
concerns the existence of the magnetic contribution to a (see section 14. 1.2[) . 

Let us for a moment deviate from approach (iia) and accept the concept employing the two quenching 
functions gx and gu- In order to get enough constraints for their unambiguous determination results of 
simulations with different (kinetic and magnetic) forcings were employed. This seems to be again ques- 
tionable, as only the structures of relations like (|64p or (|65p can be considered as universal in their validity 
ranges, but not the specific values of r* or of gx and qm- These values should in general differ for different 
realizations of turbulence. As a consequence, they can at best be determined using data from one and the 
same simulation or from simulations which are equivalent on the level of mean quantities. Exami n ing the 
visualizations of the employed realizations of turbulence given in iBrandenburg and Subramanian 



(figures 10 and 11) it is by no means clear whether this requirement is fulfilled. But even if so, an un- 
ambiguous determination of the two quenching factors were then again impossible because the relevant 
equations (17) of that paper would just become linearly dependent. 

In view of some support o f app roach (iia) from numerical studies we mention also the investigations 
of iBlackman and Field! (2002) and lField and Blackmanl (120021 ) on the saturation of a mean-field dynamo. 



Specific results on helicities obtained in this approach, even if using a relation for a like (|64p i nstea d of one 



like (|65p . are surprisingly well matched to those of numerical simulations by Brandenburgl ( 200ll ). which 
are independent of such an approach. 

Summing up, we state that a convincing numerical support of the r-approach still requires further efforts 
beyond the work published so far. 

In this context s t udies of the behavior of passiv e scal ars in a turbulent fluid carried out by 
Blackman and Field (120031 1 and by IBrandenburg et al\ l[2004l ) using the r-approach, too, deserve some 



interest. Numerical simulations confirm the results derived in the framework of this approximation rather 
well. This is remarkable since the discussed concerns with respect to the treatment of the kinetic forcing 
term (see section I4.2.2P apply here, too. 



5 An alternative 

As explained above approach (i) is based on approximative solutions of the induction equation ((4|) for b. 
The approximation is defined by the restriction to second-order correlations in the velocity fluctuations 
u. In contrast to this in the approaches (ii) neither the induction equation nor the momentum equation 
are really solved. Instead simple assumptions are used for some crucial quantities. These assumptions, 
however, have proved to be very problematic. We believe that a real progress in the calculation of the 
mean electromotive force £ beyond the second-order correlation approximation can only be gained on the 
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basis of improved solutions of the induction equation (j4]) . In the expressions for £ obtained in this way the 
quantities depending on u can then be specified according to the momentum equation. 



5.1 Higher — order correlation approximations 

As is well known it is possible to proceed from the s econd-order correlation approx imation for the mean 



As is well Known it is possible to proceed irom tne s econd-order correlation approx imation lor tne mean 
electromotive force £ to higher-order approximations (IKrause and R,adlerlll97ll . ll980h . Basic ally arbitrarily 



high orders can be reac hed. It has been shown that the proposed procedure converges ( Krausd 1968 



Krause and Radler 1980l ). However, the calculation of £ in a higher-order approximation requires the 



knowledge of higher-order correlation tensors and is moreover very tedious. That is why so far only in a 
few simple cases higher-order results have been given (e.g., Nicklaus and Stix 1988, Carvalho 1992, 1994, 
Radler et al. 1997). 

In appendix |E] a procedure is described for the calculation of and bij k in higher-order correlation 
approximations. For the sake of simplicity it is restricted to cases with U = and V • u = 0. It turns out 
that then aij and b^ k can be represented in the simple form 

/•OO !• QQ(^ f \ 

Jo Joo VQj 

b ijk (x,t) = e Um / / G(^,r)u/(x,t)w mifc (x,t;-^,-r)d 3 ^dr. (88) 

JO Joo 

The quantities v m and wij k are given as sums of constituents v^} and where v is the order in u. 

Recursions for the calculation of Vm from ^ and of from wf- k ^ are available, see appendix lEl 

Let us add a remark concerning £(°\ As explained above, equation @ for b in its original form, that 
is, without the neglect of G, with a properly specified u may have non-decaying solutions even if B 
is equal to zero. Then may be unequal to zero. As already noted in section 13.1.11 in the second- 
order correlation approximation, when ignoring the influence of the initial b, no such solutions b exist 
and therefore £^ vanishes. Hence, if we start from the second-order approximation and proceed then to 
higher approximations we never find £^°\ though it may well be non-zero. In that sense, the iteration 
procedure is not complete. 



5.2 A closure proposal 

Instead of approaches of type (ii) to the mean electromotive force £ we propose another one. For the sake 
of simplicity we restrict ourselves again to the case U = and ignore the influence of any initial value of 
b, that is, consider the limit to ~~ > — oo. Then equation (J3|) is equivalent to 

f°° f dG(£ t) ( 

b k (x,t) = (5 km 5i n - 5 kn 5im) i I — — r — kn(x - £,t - t) B n (x - £,t - t) (89) 

Jo Joo OQ v 

+(n m (x -tt-T)b n (x-tt- r))') d 3 ^ dr . 

This leads immediately to 

[•oo p dG(£ t) / 

£i(x,t) = (5ue kmn - 5 k ie imn ) / — — ({u k (x,t)u m (x-£,t-T))B n (x-€,t-T) (90) 

Jo Joo d£i V 

- (u k (x, t) u m (x - £,t - r)b n (x - £,t - r)) ) d 3 ^ dr . 
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We write for short 

£ = £ (SOCA) + £ (uub) ; (gl) 

where £( SOCA ) is the result for £ obtained in the second-order correlation approximation, that is, with 
(|90|) if there the part of the integrand with b n is ignored, and £^ uub ) the contribution to £ defined by this 
part. 

Our proposal for a new approach to £ consists in its simplest form in putting 

£ (uub) = _ f£ (SOCA) (92) 

with some factor /, which will be discussed below. Then we have 

£ = (l-f)£( SOCA K (93) 

The ansatz (I92D corresponds in a formal sense to the ansatzes (I60p or (I79D used in the approaches (ii). 
As the time r* in these cases, the factor / remains at first undetermined. In order to find rough statements 
or estimates for the factor / a look on the results of higher-order approximations reported above is useful. 
These results suggest that, e.g., in the limit q — > the factor / has the form / = c\Rm with some positive 
c\ for small Rm, and / = c\Rm + C2-Rm 2 + . . . for higher Rm. In the limit q — > oo the factor / should 
have the form / = c\St for small St, and the form / = c\St + C2St 2 + . . . for higher St. Since / — ► as 
Rm —* or St — > 0, respectively, the ansatz f)92[) trivially satisfies the requirement that higher approaches 
should include the second-order correlation approximation in overlapping parts of their respective validity 
regions. 

In this context the a-effec t calculations for the Karlsruhe experiment are of particular interest 
( Radler et al. 1997 . 19981 . 20021 ). They have been done for the case U = and steady flows u, that is q = 0, 



with arbitrary Rm. The flows have the form u = uj_ + um, where uj_ and un are, with respect to a proper 
Cartesian coordinate system (x,y,z), defined by uj_ = (u x (x, y), u y (x, y), 0) and uii = (0, 0, u z (x, y)), and 
are periodic in x and y. Two different magnetic Reynolds numbers Rm± and Rm\\ are introduced char- 
acterizing the magnitudes of uj_ and uy, respectively. If the flow described by u shows helical structures 
as, e.g., in the case of the so-called Roberts flow, an a-effect occurs such that £ = a_i_Bj_, with some 
coefficient a±, and Bj_ defined analogously to uj_; contributions to £ with derivatives of B are here ig- 
nored. In the second-order correlation approximation the result for a±, in this case called a^ OCA \ reads 
^(SOCA) _ ao j^ m ^j^ m w ith some ao depending on the geometrical properties but not on the magnitude 
of the flow. In the general case it turns out that a± = aoRm±Rm\\4>(Rm±) with the same ao as before. The 
function (j) satisfies </>(0) = 1, decays monotonically with growing Rm± and tends to z ero as Rm i — > oo ; 



further properties and numerical values of <fi are given in some of the papers mentioned (jRadler et aZ.lll997 



20021 ). When writing in the sense of ([93]) a± = (1 — f )a ( f OCA ' > we have / = 1 — cj) and therefore / — > as 



Rm± — > and — > 1 if Rm± — > oo. 



6 Summary 



The paper provides a critical view on different analytical ways to determine the mean electromotive force 
£ in mean-field electrodynamics. First some of the findings gained with the traditional approach, re- 
duced to the second-order correlation approximation or approach (i) , are summarized (section [3]) . In this 
context also some new results concerning the case of magnetohydrodynamic background turbulence are 
given (section I3.2H . Then the essentials of two versions of the r-app roach, or approach (ii), are repre - 
sented, that is, of th e simp le r- approach, or approach (iia), as u s ed by Vainshtein and Kichatinovl ( 19831 ). 
Blackman and Field! (12002') and iBrandenburg and Subramanianl (2005bl), and of t he spectral r-approach, 
or approach (iib), used bv lRadler et al. ( 20031 ) and by Rogachevskii and Kleeorin (|2003Ll2004l ) (section ED. 
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Approach (i) in its original form, applying to purely hydrodynamic background turbulence, is based on 
solutions of the induction equation for the fluctuations b simplified by some approximation. In the case 
of magnetohydro dynamic background turbulence in addition solutions of the momentum balance for u, 
again in some approximation, are used. In the approaches (ii) these equations are not really solved. Instead, 
some assumptions on crucial quantities are introduced. There is hardly any doubt in the correctness of the 
results of approach (i) in the range of its applicability, which is well defined at least by sufficient conditions. 
It is not surprising that the approaches (ii) deliver the same vectorial structures of the contributions to £ 
as approach (i), for these are already determined by elementary symmetry principles. The results of the 
two types of approaches differ however in the dependence of the coefficients of the individual contributions 
to £ on the correlation properties of the velocity fluctuations u. The discrepancies of the two approaches 
as well as strange aspects or shortcomings of the approaches (ii) are discussed in detail (sections 14.1.21 and 
I4.2.2p . Unless any overlap of the ranges of applicability of the approaches (i) and (ii) can be excluded, 
what would be very surprising, we have to conclude that the approaches (ii) are in some conflict with 
the basic equations and not all of their results can be taken for granted. For instance, the magnetic 
contribution to the a-effect occurs in the approaches (ii) not only in the case of a magnetohydrodynamic 
background turbulence but for purely hydrodynamic background turbulence, too, where they should not 
exist. In addition, in the first case it is the same "correlation time" which occurs with the kinetic and 
the magnetic contributions, what seems to be in conflict with numerical simulations (section l4.3p . Further 
the dependence of a and (3 on the magnetic diffusivity 77 and the kinematic viscosity v which occurs in 
approach (iib) is not correct. Moreover, it is doubtful whether the approaches (ii) can describe, e.g., the 
ft x J-effect correctly. 

There is no hint that the mentioned shortcomings of the approaches (ii) result from the r-approximation, 
that is, the reduction of third-order correlations of u or b to second-order ones. An important reason for 
them seems to be an improper treatment of a term connected with the forcing term in the momentum 
balance. 

A proposal is made for a new approach which avoids any conflict with approach (i) (section [5]) . In this 
context also a new formalism for the higher-order correlation approximation is presented (section I5.ip . 
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Appendix A: The limits q — > oo and q — ► 

The relations (|25p for ai v and bi pq can be written in the form 

roc „f/2 ia^i 



ejjJlp ~ e ljp 5i n n5/2 t c f f°° exp(-gg / /4rQ r , u / ,3,; ,/ 

2(4tt)3/2 Q X c LJ ^(x,aA c ,-rr c )^dedr 

eijn^ip - e ijp 5i n 3/2 t c f f 00 exp(-gg /2 /4r / ) dQ jn (x,t;(j'\ c , -t't c ) 3 , 

J r' 3 / 2 ae'i 4 



(4vr)3/2 * \ c 
2(4vr) 3 / 2 



g 5/2 r c I I ex P ( g^/^) Q j(|(x . /: g % ._/ r( , K ^ (/fl :»^ flr /. 

Joo ./0 T 



In order to evaluate (|Aip in the limit q — > oo it is useful to introduce the new integration variable 
£" = y/q£'. Then q occurs at no other places in these relations than in the third argument of Qj n , which 
takes the form ^"X c /^- Note that (,"\ c /y/Q = C"\/V T c- We may have q — > oo as a consequence of r\T c — > 
or of A c — * oo. In the first case this third argument tends to zero as q — > oo. In the second case, which 
seems rather academic, Qj n looses its dependence on this argument as q — ► oo, and we may replace the 
latter simply by zero. Considering then 

/ ^W d^(4,f, f gffl i^.ll^t, (A2) 

we find (1291) . 

In view of the limit q — ► we introduce r" = r'/g. Then g occurs only in the fourth argument of Qj n , 
that is, t" qr c = T"\ 2 c jr\. We may have q — > due to A 2 /?? — > or due to r c — ► oo. In the first case the 
fourth argument tends to zero as q — > 0. In the second one Qj n looses its dependence on this argument as 
q — ► 0, and we may replace it by zero. With 



f 

J a 



exp(- C V4r) dT= 2^ i t~ = (A3) 



r 3/2 ' C ' Jo T 5 / 2 
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we find then (1301). 



Appendix B: An alternative representation of a. and f3 as given by (|32[) 

Let us consider a homogeneous statistically steady turbulence. In this case the correlation tensor Qij 
defined by ([23}) depends no longer on x and t, that is Qij = Qij{$,,r). We introduce the quantity Uij by 



(Bl) 



POO 

Uij= / Qij(0,-T)dT. 

Jo 

With the Fourier transform Qij of Qij it holds 

Uij = J J Qij(k,-T)d 3 kdT. (B2) 

Assume now that Qij allows a separation ansatz of the form 

Q ij (k,r) = Q ij (k,0)f(k,r), (B3) 
where, of course, f(k, 0) = 1. Then we have 

Uij = J 4i(k, 0)f c (A;) d 3 k = (2tt)- 3 J Qij{i, 0) r c (0 d 3 £ , (B4) 

where f c {k) is defined by 

/>oo 

r c (k) = / f(k-r)dr, (B5) 
J o 

and t c (£) such that f c is its Fourier transform, 

Tc{0 = f Uk) exp(ik • £) d 3 fc . (B6) 
Comparing (|B 1|) and (|B4|) and remembering (|23p we arrive then at the remarkable relation 

f'OO f' 

J (n i (x,t)n i (x,t-r))dr = (27r)- 3 y ( Ui (x, t) u,(x + £, t)> r c (£) d 3 £ . (B7) 

We restrict or attention now to an incompressible fluid and assume the turbulence to be isotropic. As is 
well known the Fourier transform Qij of Qij has then the form 



Qij{k,T) = (5ij - -y L )g(k,r) +ie ijk h(k,r) . (B8) 



For the calculation of (5 according to (|3ip it is sufficient to consider a mirror-symmetric turbulence, that 
is h = 0. Then Qij satisfies the assumption (lB3j) . Consequently (1B7D applies, and (3 as given by (|32|) can 
also be written in the form 

= 37^)3 / < u ( x > *) • u ( x + *)> d3 £ • ( B9 ) 
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We may repeat the above derivations with 

Qij(€,T) = e jk i (BID) 

instead of Qij(£,r). Using again (|B8P we have then 

~ kk ■ 

Qij(k,T) = (<% - -^-)h(k,r) + ie ijk k k g(k, r) . (Bll) 

For the calculation of a it seems to be justified to put here g = 0. Then (jBll[) satisfies the analogue to 
(IB3j) , and we find an analogue to (1B7I) with Uj replaced with (V x u)j . With this result a as given by (|32|) 
can be written in the form 

« = / < u (x, *) • ( V x u ( x + *> *))> Tb(0 d 3 e • (B12) 

The quantities r c in (|B9j) and (|B12h may be well different from each other. 



Appendix C: tt X J effect 

Under the assumptions which justify equation ([8|) we have 

hjk = ^Cijfc ~~ <^l% ^k — ^2^ik % — S^Sjk f2j . (CI) 

The /3, <5i and (52 terms lead just to the last line of ([8]) whereas the 83 term has because of V • B = no 
counterpart there. Multiplying both sides of ()Cip first by 8ij fi^ and then by <5j/% Clj and by fij we arrive 
at a system of three equations for <5i, 82 and £3. Its solutions for 8\ and £2 read 

1 1 
£1 = - 1Q ^ 2 (4&iij - biji - bju) VLj , S 2 = {buj - Abiji + bju)Qj . (C2) 

Expressing 6^ according to ([25]) we see that 81 depends only on the part of Qij which is firstly an- 
tisymmetric in the indices and secondly even in Assuming that the turbulence is homogeneous and 
statistically steady, which implies Qij(£, r) = Qji{—£, — r), we may conclude that this part is odd in r and 
therefore vanishes at r = 0. As for 82 this statement applies only under the additional assumption of an 
incompressible fluid, that is, dQij/d^j = 0. 

Thes e conc lusions can also be drawn immediately from the first results on the fi x J-effect given by 
Radlerl ([1969). By the reason explained above, the functions h and k used there are odd but I is even 



in r, and it can be easily shown that / vanishes in the incompressible case. This just confirms the above 
statements on the part of Qij responsible for the coefficients that correspond to 5± and 82- 

By the way, the rea son why the ft x J-effec t did not occur in the first calculation of the mean electromotive 
force £ delivered by ISteenbeck et al. was a not completely correct ansatz for Qij, just without a 

part being odd in r. 



Appendix D: Derivation of equation (|68p 

We start from (|67p and express vt and Wj on the right-hand side by their Fourier representations, 

<Mv,w;x,£)= ! [{viWwjOs.")) exp(i(k'-(x + ^/2)+k"-(x-^/2)))d 3 A ; / d 3 A : ,/ . (Dl) 
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Introducing new integration variables k and K by k' = k + K/2 and k" = — k + K/2 we find then 



*ij(v,w;x,£) = J 6 ij (x,k)exp(ik-^)d 3 A: 

^■(v,w;x,k) = J(v i (h + K/2)w j (-k + K/2)) exp(iK • x) d 3 ^ . (D2) 



The last line is identical to (1681). 



Appendix E: Higher order correlation approximations 

Consider again an infinitely extended homogeneous fluid and assume for the sake of simplicity U = and 
V • u = 0. Assume in view of the calculation of £ at x = x that 

Bi = Bi + Bijixj-Xj), (El) 

with Bi and By independent of position and steady, and Ba = 0. Clearly Bi and Bij are equal to Bi and 
dBi/dxj, respectively, at x = x. From and ([5]) it follows that 

d t b k - r]Ab k = (B p + B pq (x q - x q )) - B kp u p + (5 kq 5 pr - S kr 5 pq )-^-(u q b r )' , = (E2) 

and, if again the influence of the initial b is ignored, that is to ~~ * — oo, 

b k (x,t)=B p HI U k ( X -tt-T)d 3 tidT 
JO Joo C^p 

" B pq ri (g<"> (£, r) 5 kp 6 gr - dG{n) f^ ( x _ £ _ f ) «5 fcr ) Uj .( x _ £ t _ T ) d 3 £dr (E3) 

f°° f dG^) (£ r) z 

Jo Joo 

Put then 

h = b[ l) + bf +bf + ... (E4) 

and obtain 

Jo Joo 

-^w / °°/ (£, t) 6 kp S qr - 9G " & T) (x q -x q - e a ) flJ u r (x-S,t- t) d 3 £dr (E5) 

J0 Joo V 7 

6^(x,t) = (5 kq 5 pr - 5 kr S pq ) HI' dGiV J C,T \ u q (x -tt-r) b^\x -£,t- r))'d 3 £dr 

Jo Joo °V> 



v>2. 
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These relations are equivalent to 

JO Joe a ?p 

+B pq f°°f G^\^T)(w^(^t-,-^-T))'d^dr, u>l, 

JO Joo 

with 



4 1} ( x '^^' r ) = n fc(x + £,t + t) , 

wJb (x,*;€,r) = / / —7 U fc (x-£ ,t-t)^ (x - £ , i - r ; £, r) 

JO Joo v 



and 



-u, (x - r - r') ] (x - * " t';£, t) d^'dr' , 



W fepg( X '*;€' T ) = -4p«q(x + |,t + r) + (X q - X q + £ g ) - -> ' ' 



3£ P 

«£° (x, t; r) = jTj^ ^^f^ («*(* - r - r') ^ (x - t - r'; r) 



fcpg 

From (|E4|) and (|E6[) it can be concluded that 



Ul (x - r - r') (x - t - r'; r) ) d 3 £'dr' . 



b k (^t)=B p [°°[ dG l^ T \ v k ^,t;-t-r)yd^dr 

Jo Joo C?p 



P 

oo 



+B pq / / GW^rJC^x.tj-e.-rj/d^dr, 

JO Joo 



where 



w fc (x, t; £, r) = 4 1} (x, t; £, r) + (x, t; £, r) + 4 3) (*, *; £ , r) + . . . 
w fepg (x, t; t) = wjjj^x, t; £, r) + wj^x, t; r) + wjj ? (x, t; £, r) + . 

with and as given by (|E7p and (|E8|) . 

Calculating then Si = €ijk{ujbk) and expanding it according to © we finally find 

— - — Uj{% t) u fc (x, £; -t) d 3 £ dr 

_o C?p 

6 ip5 (x, t) = e ijk / / G M (£, r) «j (x, t) w kpq (x, t; -r) d 3 £ dr , 

JO Joo 

which has to be completed by ([E7j) . fjE8|) and (jE10|) . 
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Note added in proof 

In appendix [B] it has been shown that the results (|32j) of approach (i) for a and 0, with integrals of 
(u(x, t) • (V x u(x, t — t))) and (u(x, t) • u(x, t — r)) over r, can be rewritten into (1B9|) and (|B12j) . with 
integrals containing u only as (u(x, t) ■ (V X u(x + t))) or (u(x, t) ■ u(x + t)), taken over all £. We have 
argued there with very specific properties of the turbulence, that is of u. 

Relations between the two types of such integrals, over r and over exist also in a much wider range 
of assumptions. Consider any function F = F(£, r). Assume that the integral F(0, — r) dr has a non- 
zero finite value. Assume further that the integral J F(£, 0) g(£) d 3 £, with any weight function g, has a 
non-zero finite value, too. Clearly the magnitude of a given g can be fixed such that 

co r 

F(0,-r)dr = (2vr)- 3 / 0) d 3 £ . (*) 

J CO 

There is, of course, some freedom in the choice of the dependence of g on We may determine a specific 
g, e.g., after the pattern of appendix [Bl For this purpose we introduce the Fourier transform F of F with 
respect to £ and put F(k, — r) dr = F(k, 0) h(— k). Integrating both sides of this relation over all k 
and using the convolution theorem we arrive just at (★), with g(£) such that its Fourier transform is equal 
to h(k). 

Let us first identify F in relation (*) with (u(x, t) ■ (V X u(x + £, t — r))) or (u(x, t) -u(x+£, t— t)}. Then 
we arrive without utilizing specific properties of the turbulence at relations which allow us to rewrite the 
results ([32} for a and f3 into §B9§ and (|BT2|) . In the same way we may rewrite the more general results 
(|29|) for dip and bi pq such that, as in (f85|) . only correlations of components of u at the same time and only 
integrals over £ occur. 

This finding extenuates our statement made in section 14.2.21 on the discrepancy between the results (J29]) 
of approach (i) and the results (|85p of approach (iib). However, it probably not even completely resolves 
the part of the discrepancy resulting from the different types of correlations and different types of integrals: 
So it is hard to imagine that all integrals in (|29[) over correlations of u can be rewritten as integrals over 
£ with one and the same g, which would have then to be identified with the r* that occurs in (I85p . 
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